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These notes are a summary of our efforts to answer the question in the title. Our answer is 
in the affirmative as: (i) HERA data indicate a large value of the gluon structure function; (ii) 
no contradictions with the asymptotic predictions of high density QCD have been observed; 
and (iii) the numerical estimates of our model give a natural description of the size of 
deviation from the routine DGLAP explanation. We discuss the alternative approaches and 
possible new experiments. 



I. PROBLEMS: 

In the region of low x and low Q 2 which is now investigated by HERA we face two challenging problems : 

1. Matching of "hard" processes, that can be successfully described using perturbative QCD (pQCD), 
and "soft" processes, that should be described using non-perturbative QCD (npQCD); 

2. Theoretical description of high density QCD (hdQCD). In this kinematic region we expect that the 
typical distances will be small but the parton density will be so large that a new non perturbative 
approach needs to be developed for dealing with this system. 

II. MAIN IDEA: 

The main physical idea, on which our phenomenological approach is based is JiJ : 

The above two problems are correlated and the system of partons always passes through 
the stage of hdQCD ( at shorter distances ) before it proceeds to the black box limit, which 
we call non-perturbative QCD and which, in practice, we describe using old fashion Reggeon 
phenomenology. 



III. SCALES IN DIS: 



Since HERA started to investigate a new kinematic domain of hdQCD, we want to find the most funda- 
mental phenomena which are typical of hdQCD as well as npQCD in the HERA experimental data. We 
would first like to determine what the values of the scales or typical distances at which two transitions 



pQCD — ► hdQCD — > npQCD 

occur. However, before discussing these transitions let us recall that a photon - hadron interaction at high 
energy in the rest frame of the target has two sequential stages in time:(i) 7* — > hadron system ( qq - pair 
); and then (ii) this hadron system ( qq - pair ) interacts with the target |2|. This time sequence allows us 
to write the cross section for photon-proton interaction in the form: 

atot{l*p) = ^2 l^nf <rtot(n,x) , (1) 

n 

where ^ n is the wave function of the hadron ( parton ) system produced in the first stage of the process. 

The first scale we introduce to separate the long and short distances is r s ^ p . Roughly speaking, the QCD 
running coupling constant as(r^) is small as(r^) < 1 for r± < r s ^ p , while as{r\) rj 1 for r± > r± p . 

The second scale we associate with the transition between the low parton density phase and the high parton 
density phase of pQCD. This scale depends on the energy of the colliding system and can be determined 
from the condition that the packing factor (PF) of partons in a parton cascade is equal to unity |l| : 

_ 3 7r 2 a 5 xG(x,Q 2 s (x)) 
PF = K = oTTSTZT x = !; ( 2 ) 



2Q 2 s (x) irR 2 



r saturation ~ ^/Qs( X )- 



Fig. 1 shows this packing factor in the HERA kinematic region. One can see when both Q and x are 
sufficiently small, we really have a dense parton system at HERA. 




FIG. 1. The packing factor for the gluon from HERA experimental data. 



IV. THEORETICAL STATUS OF HIGH PARTON DENSITY QCD: 



In DIS at low x one can find a high density system of partons, which is a non-perturbative system due to 
high density of partons although the running QCD coupling constant is still small ( as{r±) « 1 ). Such a 
unique system can be treated theoretically |J . It should be stressed that the theory of hdQCD is now in 
very good shape. 
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Two approaches have been developed for hdQCD. The first one (|] is based on pQCD ( see GLR and 
Mueller and Qiu papers in Ref. jjj ) and on dipole degrees of freedom [Q. This approach gives a natural 
description of the parton cascade in the kinematic region for k < 1 and up to the transition region with 
K « 1 ( see Fig. |). 



ln(l/x) 

A 



K» 




K = 1 



K « 1 



000 



FIG. 2. The parton distribution in the transverse plane. The curve shows the saturation scale Q s (x;A) 

The second method uses the effective Lagrangian suggested by McLerran and Venugopalan |j| , this is a 
natural framework to describe data in the deep saturation region where k 3> 1 (see Fig. |^) . As a result of 
intensive work using these two approaches the non-linear evolution equation has been derived Q which has 
the following form 



da el (x 01 ,bt,y) 
dy 



2Cf a s 



In 



a el (x,b t ,y) + 



C F a s 



x 2 



2 a' 



; (x 2,bt - -Xi2,y) - 



'(x 2,b t - -Xi2,y) a e '(xi2,b t - -x 2,2/) 



(3) 



where a el (r^, b t , x) is the elastic scattering amplitude of the dipole of size r± at energy ocl/x and at impact 
parameter bt- The cross section in Eq. ([I]) is equal to a(r^,x) = 2 J d 2 bta el (r^,b t ,x). The pictorial form 
of Eq. (^) is given in Fig. ^ which shows that this equation has a very simple physical meaning: the dipole 
with size x\o decays in two dipoles with sizes X\2 and 202 ■ These two dipoles interact with the target. The 
non-linear term which takes into account the Glauber corrections for such an interaction, Eq. (||) is the 
same as the GLR -equation jlj but in the space representation. It gives a correct coefficient in front of the 
non-linear term which coincides with one calculated in Ref. & in the double log limit. We wish to stress that 
this equation which includes the Glauber rescatterings , has definite initial conditions and has been derived 
by both methods (see Refs. |||8|]). 

We have devoted much time to the discussion of the pure theoretical approach as a comparison of our 
model approach with the solution to Eq. (g) |7],|9j which is the criterion of how well or how badly our model 
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works. However, we still need a model because (i) we have to choose a correct initial distribution to solve 
Eq. (|J); (ii) to determine the value of the phenomenological parameters which enter Eq. (J3J) through the 
initial conditions; and (iii) to find out and extract the value of the gluon structure function from a{r\, bt,x). 
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FIG. 3. The pictorial form of the non-linear evolution equation in the hdQCD kinematic region. 



V. MODEL: 



The first thing we have to specify when constructing a model, is the choice of the correct degrees of freedom 
(DOF) or, in other words, we have to answer the question what we mean by letter n in Eq. ([j]). 

It is well known that for short distances (r± < r s ± p ) the correct degrees of freedom in QCD are colorless 
dipoles Q] . This means that in Eq. ([j]) at short distances 

*n — > *(Q 2 ;rl,z) ; (4) 
a(n,x) — ► a(rj_,x); (5) 

where Q 2 is the photon vituality, z is the fraction of energy carried by quark and r± is the dipole size Q 

However, at long distances (rj_ > r s ^ p ) what the correct degrees of freedom are, is still an open question. 
In our approach we use the constituent quarks as the degrees of freedom. This assumption leads to 

*n > ^hadronin) \ (6) 

a(n,x) — ► cr q (xi); (7) 

where Ti are the quark coordinates and a q {xi) is the cross section of the constituent quark i with fraction 
of energy Xi with the target. This assumption is certainly a pure model assumption and an argument in 



1 Actually, Eq. (Q) and Eq. || were understood long before Ref. Q. In Ref. |ll[ it was shown that the dipole size is 
the correct degree of freedom for the interaction of quark-anti quark pair at high energy, if this interaction is induced 
by two gluon exchange. In Ref. jl2| this claim was generalized for the leading log approximation of pQCD. In Ref. 
[ p^[ it was proved that a dipole picture can be used for describing the gluon - hadron interaction. However, the final 
point in Ref. [^] was the most important since it was proven that the QCD evolution at low x can be rewritten in 
terms of colorless dipoles. 
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support of it is the fact that the same assumption is a part of the successful Regge jl4| and VDM jla 
phenomenology and that such degrees of freedom appear in the instanton models of the QCD vacuum Mm 
Therefore, for long distances the scattering amplitude Eq. (Q) can be written 



Vtot{l*P) 



E 

M„ <Mn 



Q 2 + M2 



i*nr<)i a E ff *fo) ' 



(8) 



where Mq w l/r s ^ p . It is clear that Eq. (||) is nothing more than the VDM model with AQM prescription 
for the rescatterings of vector mesons with the target ( see more details in Ref. JEtJ ). 
For short distances ( r± < r^ p ) Eq. (0) reduces to the form 



Vtot(l*P) = I d 2 r ± / dz |* 7 {Q 2 ;r ± ,z)\ 2 <j d ipoie(r±,x) , 



where the wave function of the virtual photon are well known |13| , |18| . 
For (Tdipoie ( r ± ' x) we use the Glauber - Mueller (eikonal ) formula |y] [ll 



where 



■ole(r±,x) = 2 



Cl(r±,x; b t ) 



d 2 b t { 1 



4 

xG(x, —;b t ) 



(9) 



(10) 



(11) 



3nR 2 ■ 

It should be stressed that in our calculations we used for G(x, -4-;M the Glauber-Mueller formula ||l3|| 

r ± — 

namely, 

AttR 2 f 1 daS_ r°° dr' 2 

x 



with 



4 

xG(x, —;b t ) 



r' 4 



n 



DGLAP 



(12) 



(13) 



Table 1 


Reaction 


Q 2 (GeV 2 ) 


X 


References 


(Ttot{-y*p) 


4-65 


< 0.01 





F2(x,Q 2 ) 


14-65 


< 0.01 




xG{Q 2 ,x) 


14-65 


< 0.01 


§ 


dF 2 /dlnQ 2 


14-65 


< 0.01 


urn 


fft ot (77*) 


0;0 4-20 


< 0.01 


n 


diff 
°tot 


5 4-65 


< 0.01 


11 


diff 
"tot 
"tot 


14-65 


< 0.01 


11 


cr(7*p^ J/^+p) 


4-65 


< 0.01 


@§ 


slope B(-y*p — > J/$ + p) 


4-65 


< 0.01 


@§ 


slope B(-y*p — > p + p) 


5 4-65 


< 0.01 


§3 
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In Eq. (|13|) we use the Gaussian parameterization for the profile function for the quarks in the proton 

where R is the proton radius. 

Eq. ([h]) - Eq. (|l4|) take into account the rescattering quark - anti quark pair and one ( the fastest) gluon in 
the parton cascade. It turns out that the rescattering of the fastest gluon is very important in our attempts 
to describe the experimental data. The data, which we described, as well as the details of our model, can 
be found in references, presented in Table 1. 

The main features of our model are summarized in the following Table2. 



Table 2 



Perturbative QCD 
short distances 

r±_ < 

DOF: color dipoles j| 

ty„: QED for virtual photon 

Gtotin, x) — a(rt,x) 

Mueller-Glauber Eikonal O for a(rl,x) 



_L 



non-perturbative QCD 
long distances 

< r± 

DOF: constituent quarks ^) 

ty-n'- generalized VDM for hadronic system 

o t ot{n, x) — a(qq — > qq; x) 

Regge phenomenology for a(q + q —> q + q;x) 



Fig. U shows the r± dependence of the total and diffractive dipole cross sections, where crfof ole is calculated 
using Eq. ( |Io| ) - Eq . (|l4| ) and where the diffractive cross section for a dipole is estimated by the Kovchegov- 
McLerran formula 129(1 
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FIG. 4. The r± dependence of uf^ ole (Fig. 4~ a ) an d ct^? ° f te (Fig.^-b). 
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VI. ADVANTAGES OF THE MODEL: 



For short distances r± < r satura tion ~ ^/Qs{x) the model agrees with the pQCD predictions. They 
are the same as in DGLAP evolution equations; 

The model gives a natural explanation and an estimate for the value of the saturation scale (see Fig. |5|); 




X 



FIG. 5. The saturation scale Q s (x) in our model . 

This model is a good approximation to the full solution of the nonlinear evolution equations that 
describe theoretically the hdQCD kinematic region ( GLR -equation and its' generalizations |l|||^]|6)). 
See Figs.6-a - 6-c and discussion below; 

Our model correctly reproduces the operator product expansion (OPE) and gives the higher twist 
contributions with the correct anomalous dimension |p7f in the region of low x; 

In the model as well as in Mueller- Glauber approach in general one preserves the relation between 
elastic and quasi-elastic scattering and multi particle production in DIS based on the AGK cutting 
rules p6[ . This allows us to describe the inclusive cross sections as well as the correlation functions; 

Our model includes the impact parameter behavior of the dipole-target amplitude, consequently it can 
be generalized to DIS with nuclei. We have completed some calculations for DIS with a nuclear target 
p8| but will not discuss them here; 

The model violates the energy sum rules but the discrepancy is very small as can be seen in Fig. 6-d; 
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k values 



k values 





Fig. 6-c Fig. 6-d 

FIG. 6. Figs.6-a - 6-c show the packing factor k for gluons in the DGLAP evolution equation (DGLAP), in 
Mueller- Glauber Eikonal approach (MODMF), in the correct non-linear evolution equation (AGL) and in the simplified 
version of the GLR equation (GLR). Fig. 6-d displays the energy-momentum sum rules in our model. 



• In the region of large r± > r s ^ p the model displays all the attractive features of the VDM and is able 
to describe "soft" physics in DIS well. 



VII. DISADVANTAGES OF THE MODEL: 



We have two major problems with this kind of a model. First, the theoretical accuracy of the general 
formula 

<Jtot(j*p) = J d 2 r ± dz \^"(Q 2 ;r ± ,z)\ 2 a dipo i e {r 2 ± ,x) (15) 

is only approximate. This formula can be proven |I^,[l3|,^| only in the leading log(l/x) approximation of 
pQCD in which we consider as m(l/x) w 1 while as <C 1. Therefore, Eq. (||) is much less general than the 
DGLAP approach and our results are always worse than those for the DGLAP evolution equations. This is 
the price we pay for a far more transparent formalism which includes the long distance physics. 

Second, the separation scale r^_ ep changes for different reactions. Even for longitudinal and transverse 
polarized photon we have different r^ p ss 1/(0.85 ~ 0.95)GeV while r s ^ p > 0.6GeV. This fact limits the use 
of the attractive factorization properties of Eq. (|l|) or / and Eq. (||) . 

Other shortcomings: 
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• In the region of low Q 2 the model's results are always less than those found by solving the non-linear 
evolution equation (see Fig.6-a). Therefore, the estimates obtained from the model are below the 
measured effects. Knowing these limitations we plan to improve our model; 



• The eikonal approach cannot be correct at very high energies so our model is specially constructed and 
suited for the HERA kinematic region and can only be used with a lot of caution to higher energies 
and, in particular, to the LHC energy range; 

• We used a Gaussian profile function in impact parameter (b t ) which cannot describe the region of very 
large momentum transfer. Therefore, to be accurate we have to try alternative functions for large t 
processes. 

VIII. ASYMPTOTIC PREDICTIONS: 

In our model the dipole cross section at high energy and fixed r± is equal to 

o dlpole (r\,x) = 2TT<b 2 t {x)> (16) 

where < b 2 > is the average impact parameter in dipole-target scattering . As function of energy it behaves 
in our model as 



< b 2 > = R 2 (C + HQ 2 (x)r 2 ± }) 



o 2a' e f f ln(l/x) 



where R 2 is the target size and C is the Euler constant. Fig. 7-a and Fig.7-b show that the effective slope 
is rather small ( it is at least four times smaller than the soft slope of Pomeron trajectory). However, for 
the processes induced by the gluon structure function this effective shrinkage could be visible and it has to 
be taken into account. 



<b, > F , 



• Q'= 2.5 GeV 

= 5.0 GeV 2 
= 10. GeV 2 

• = 20. GeV 2 



10.0 15.0 
y=ln(1/x) 



<b, >,r 





R 2 = 10 GeV' 2 














sssl ^^^~-~ 










Q'= 2.5 GeV 






= 5.0 GeV 2 






=10. GeV 2 






=20. GeV 2 











10.0 

y=ln(1/x) 



Fig. 7-a Fig. 7-b 

FIG. 7. Figs. 7-a - 7-b show the effective shrinkage of the diffraction peak in our model for F2 and xG 



Such a behavior of the dipole cross section results in the following predictions for the observables: 



• F 2 



^YH S Z)Q 2 - <b 2 >^ Z 2 Q 2 -[R 2 + 2a' e// ln(l/x)j; 



dF 2 
dlnQ 2 



R 2 (HQ 2 (x)/Q 2 )) 



\n(Q 2 Jx)/Q 2 ) , 



2 l+ln(Q2(a;)/Q2) , 



The ratio of the diffractive to the total cross section should not depend on energy 



rdiff 



a tot 



Const 



slowly 1 



• The energy behaviour of the diffractive cross section is determined by short distances r± rj 

saturation = 1/Q g ( x ) ; 

• The high density QCD effect should be stronger in the diffractive channels. Therefore, we expect to 
see these effects firstly in diffractive production of heavy mesons or in inclusive diffractive production; 

• We expect minima in t distribution of the diffractive production (especially in the case of diffractive 
production of mesons in DIS ) which is related to wave propagating picture described by eikonal 
approach. 

IX. PHENOMENOLOGICAL PARAMETERS OF THE MODEL: 

Before discussing the application to HERA data we list the parameters that we use to fit the data. 

A. R 2 - size of the target. 

The size of the target enters the impact parameter profile of the target which we take in the Gaussian 
form: 

S(fc) = ^e-S. (17) 

The HERA data for photo production of J/ 'J - meson as well as CDF data on double parton cross section 
leads to the value of R 2 = 5 -r 10 GeV- 2 . We use R 2 = 5 GeV~ 2 and R 2 = 10 GeV- 2 for an estimate of 
the possible effect and R 2 as a fitting parameter for the description of the experimental data. Note, that the 
value of R 2 — 8.5 GeV~ 2 was taken for all reactions that we have described. 

B. Qo — l/ r se P - separation parameter. 

As we have discussed we can trust our model for the saturation effect ( see Eq. ( [U)| ) - Eq. ( 53 ) ) only at 
rather small distances (r± < r s ^ p ) or , in other words, at large virtualities of the incoming photon Q 2 > Qq. 
We have commented on the value of r^ ep , but in practice we used Qq = 0.6 -r- 1 GeV 2 and tried to study 
how our fit depends on the value of Qq. Therefore, the result of our calculations should be read correctly, 
as "the shadowing corrections from short distances r± < I/Qq gives this or that 

C. Solution of the DGLAP evolution equations. 

We tried to use all available parameterization of the solution of the DGLAP evolution equations |H]|3l| , but 
we prefer the GRV parameterization p3] . The reason for this is very simple: the theoretical formulae, that 
are the basis of our model, were derived in double log approximation of pQCD and the GRV parameterization 
is the closest one to the DLA. 

X. OUR MODEL VERSUS HERA DATA: 

A. <rtot(7* + p) at low x and Q 2 . 

As has been discussed the matching between "soft" and "hard" processes is not very sensitive to the 
saturation scale since we will show in the next section that the shadowing corrections to F% are rather small 
0. Therefore, we use a to t{l*p) to extract the separation scale from the experimental data j3^j34|]. It turns 
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out that 
0.6 GeV~ 



sep 



1/Qot ~ 1/(0.85 ~ 0.95) GeV for the transverse polarized photon, while r ± p « 1/Qql > 
Fig. H shows our description of the experimental data at x < 0.01 




1 10 

Q z (GcV 2 ) 



FIG. 8. <Tt ot (7* +p) m our model. From these data we extracted Q\ T = 0.7 — 0.9 GeV 2 and Q\ L < 0.4 GeV 2 . The 
vertical line corresponds to x = 0.01. 



B. F 2 . 



Our main conclusion that the global features of F2 is not very sensitive to the contributions of the SC 
from distances r± < r^ p 0. In Fig. || we plot the result of our estimates for low Q 2 — 1.5 -j- 6.5 GeV 2 , for 
higher Q 2 the shadowing corrections are even smaller. 
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X 

FIG. 9. F2 with ( dashed line) and without SC in our model. SC has been taken into account for distances shorter 
than the separation scale (r± < r s ^ p — 1/Qo). 

C. xG(Q 2 ,x) 

Our calculations for xG{x,Q 2 ) has been discussed in section VI ( see Figs. 6-a - 6-c and Ref. J7J ). One 
can see from these figures that unlike the case of F2 , the shadowing corrections are rather large and they 
tame the increase of the gluon density. Actually, all indications of high density QCD effects, that we see in 
the HERA data and will discuss below, are due to strong screening in the gluon channel. 



D. F 2 Slope. 

We consider the Q 2 behavior of the F2 slope, d f^Qi , as the best experimental indication of the strong SC 
or other effects of high density QCD. Our study (2Cjj25|] shows that none of available parameterization of 
the DGLAP equation can reproduce the experimental data, while our model does this. More than that, two 
different parameterization such as GRV'94 and GRV'98 both describe the data well after taking into account 
SC in the framework of our model. Fig. [To] shows the comparison of the GRV'98 parameterization with the 
HI data. One can see that GRV'98 which was invented to describe ZEUS data ( Caldwell-plot) failed to fit 
HI results but in our model we can fit HI data without any additional change of parameters. 
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Fig. [Tl] presents the 'ideal' Caldwell plot from the point of the gluon saturation (see section VII): the 
F2 slope as function of Q 2 at fixed x. One can see two features in this figure: (i) our model describes the 
HERA data quite well; and (ii) the data as well as our calculation do not show any maximum which could 
be interpret as a saturation scale. Our model is also able to describe the Q 2 behaviour of the F 2 at fixed W 
which shows maxima. It means that the maxima in the fixed W plot does not reflect the saturation phase 
of the parton cascade but are an artifact of the particular choice of the kinematic variable. 




FIG. 10. F-i-alope in GRV'98 and in our model with GRV'98 input for gluon structure function. Data are taken 
from figures of Ref. Mw and from Ref. J^ /. 
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FIG. 11. The Fi slope at fixed x versus Q .Data are taken from figures of Ref. \3by and from Ref. 



We have also studied an alternative approach, namely, the Donnachie-Landshoff model for matching the 
"soft" and the "hard" interactions J37J: the sum of "soft" and "hard" Pomerons (see Fig. |lj). One can see 
that this model is as successful in describing the data as is our model. This means that at the moment we 
cannot claim that the data show a saturation effect. We have a much more modest claim: The current data 
can be described in two different approaches, either due to gluon saturation or due to matching between "soft" 
and "hard" contributions ( "soft" and "hard" Pomerons ) but at rather large momenta ( about 1 — 2 GeV ). 
In spite of the fact that we failed to arrive at a unique conclusion one can see that measuring the F% slope 
gives new and interesting information about QCD dynamics. 



14 




ZEUS Preliminary 
10- 3 




0.1 


1 


10 


Q 


2 


0.1 


1 10 Q 2 


X = 




10" 3 


10" 2 


10 


-1 


X = 


= 10- 3 


10" 2 10" 1 


0.6 


I I I 


I 


I 1 1 

W=30 


\ 




0.6 


JHI| I 


1 1 1 1 '- 

W=14 - 


0.4 












0.4 






0.2 






***** 






0.2 













1 . 











0.1 


1 


10 


Q 


2 


0.1 


1 10 Q 2 



FIG. 12. The F2 slope at fixed W versus Q . Curves are calculated In Donnachie-Landshoff model IfiTj . Data are 
taken from figures of Ref. pal and from Ref. p4 



E. Energy dependence of the inclusive diffraction cross section cr*f and the ratio . 



In a saturation approach the typical distances that are dominant in the diffractive production, are of the 
order of the saturation scale r s ^ turatlon |l| ; ^ 39| , unlike the total cross section where they are of the order of 
(l/<5).n One can see this directly from FigT^ which shows that in the saturation region a hadron appears as 
a diffraction grid of size r s ^ turatlon , This fact is in a perfect agreement with the HERA experimental data 



58[ as well as with our estimates 22 



However, we cannot reproduce in our model ]23| ] the experimentally observed fact that the ratio a^J^ ja to t 
is almost energy independent. The dependence on energy in our calculations stems partly from distances 
shorter than 0.3/m (see Ref. |2^] for details ) where we can trust the DGLAP evolution equations. The 
DGLAP evolution was included in our model in contrast to the Golec-Biernat and Wusthoff model |4t| which 
includes the saturation scale but has a very oversimplified behaviour at short r± {udipoie ^ r ± )■ 
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FIG. 13. a dt ^ I sigmatot vesus energy in our model. Solid line is total contribution, dashed line is qq and dotted 
line is qq + G contributions. 



F. Energy dependence of J/\E' production. 

We claim that the energy behaviour of the exclusive vector meson production, especially, J/'F-production 
is very discriminating with respect to SC (see Fig. [l4|) . In spite of large errors related to uncertainties due 
to our poor knowledge of the wave function of vector mesons, this uncertainty contributes mostly to the 
normalization of the cross section while the energy slope is still a source of the information on the SC. We 
believe that at the moment there is no parameterization of the DGLAP evolution equation which is able to 
describe simalteneously the i^-slope and the energy dependence of the J/'F photo production where a huge 
amount of new data has been accumulated ( see Ref. p5| and references therein for more details ). 
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200 




FIG. 14. Examples of the description of J/9 DIS and photo production. We got excellent n.d.f. < 1. 

G. t - behaviour of hard diffractive leptoproduction of vector meson. 

As has been mentioned, in our model we generate a shrinkage of the diffraction peak (see Fig. ^ ) as well 
as a damping of the energy dependence. We want to draw the readers attention to the fact that we expect 
a characteristic t-dependence with possible minimum around \t\ = 0.8 -j- 1.2 GeV 2 . 
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fig.l5-a Fig. 15-b 

FIG. 15. The energy behaviour of the diffraction slope ( Fig. 15-a) and the t - dependence (Fi g. 15-b) for J/ty 
photo production. Fig. 15-b corresponds to the two radii model for the impact parameter dependence \2J^ of the profile 
function which provides a better description of the experimental data. Fig. 15-a is plotted for one radius model. The 
experimental data are taken from J^_?|/. 



XI. SEARCHING FOR NEW OBSERVABLES: 



A. <x(7 7*). 



The interaction of two virtual photons provides a unique opportunity to measure both the property of the 
BFKL Pomeron [B as well as to explore the saturation region. By now we have checked that our model is 
able to describe plf current experimental data for the jj* total cross section [|43| with the same parameters 



that we used for DIS with the nucleon target (see Fig. 16 ). However, we cannot reproduce in our approach 
the experimental data for the cross section of two virtual photons with Q\ k, Q\ > 1 GeV 2 . Perhaps, we 
need to improve our approach including the BFKL contribution. 



B. Maxima in ratios. 



In our attempt to find an improved observable which will be more sensitive to the saturation scale we 
study the Q 2 behaviour of the ratios: Fl/Ft and F® /F® for longitudinal and transverse structure function 
for inclusive DIS and for diffraction in DIS 0]. We found that these ratios have maxima at Q 2 — Q^axix) 
which are moving as a function of x. In Fig. ^ we have plotted some examples of these ratios and the 
behaviour of Q max (x). It appears that Q m ax(x) is a simple function of the saturation scale Q 2 (x) of Eq. (^). 
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FIG. 16. The description of 7* 7* total cross section at different values of the photon vituality in our model. The 
experimental data are taken from Ref. j^g/ . 




C. Higher twist contributions. 



As we have mentioned, one of the attractive features of our model is the fact that we are able to describe 
the higher twist contribution in accordance with all known theoretical information. Our calculations confirm 
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the result of Ref. |4jj] that there is an almost full cancellation of higher twist contribution in F 2 while they 
give substantial contributions separately in Fl and Ft as well as in F D . We hope, that the twist analysis 
will help us to find a scale when the higher twist contribution will be of the order of the leading one. It is 
well known jjj, that in a saturation approach at Q 2 — Q 2 (x) all twist contributions are of the same order. 
Therefore, our twist analysis will assist in finding a kinematic range of Q 2 and x where the higher twist 
should be important and where we expect a manifestation of the gluon saturation. 



XII. OUR PLANS FOR THE NEAR FUTURE. 

We list here the problems that we have started to work on: 

1. Large pt-jet production in photo production. The first calculation indicates a maximum at the satu- 
ration scale in p t -spectrum; 

2. Heavy quark production and especially open charm production will appear soon; 

3. Using our model we intend to calculate the inclusive observables and to examine their sensitivity to 
SC; 

4. We have started, and we are going to continue, to check our approach with the 7*7*- data on total 
cross section; 

5. The calculations for gluon structure function as well as for F2 for nuclear targets has been performed, 
as well as predictions for the total diffractive cross sections for DIS with nuclei. We have extended 
our twist analysis to DIS with nucleus |2^] to see whether we will be able to separate the higher twist 
contribution using A-dependence. We have started to calculate all above mentioned observables for 
DIS with a nuclear target, hoping that the future experiment for DIS with nuclei will open a new 
dimension in our investigation of the saturation phenomenon and will be complementary to HERA 
data. 



XIII. FEEDBACK FOR THEORY: 

It should be mentioned that our model or a similar type of approach is needed to solve the nonlinear 
evolution equation that have been discussed in section IV. The essential requirements to solve Eq. (||) are: 

• a correct initial distribution, which should be of the eikonal type as was shown by Kovchegov Q; 

• values of the parameters that enter Eq. (|^), such as the 6 t -distribution and a separation scale; 

• a relationship between the dipole scattering amplitude a el {r\, bt', x) and the gluon structure function. 

Our model provides an answer to all formulated questions and can be used as an initial distribution for 
the nonlinear evolution equation (see Eq. ). Indeed, our model gives the eikonal description of the 
experimental data which is suited to be used as the initial condition for Eq. (|^). The only item that we 
need to add is the relation between the gluon structure function and the dipole scattering amplitude. This 
relation is a natural generalization of Eq. ( |l2|) , namely, 

xG(x,—) = / d 2 b t / — / -A 2 Im a{rl,b t ;x) . (18) 

r ±_ I? J Jx x Jr ± r ±_ 

A similar relation for the F2 structure function stems directly from the relationship between the measured 
cross section ( a(j*p) ) and F 2 ( cr(~f*p) = ^p-F 2 ) and from Eq. (||) and Eq. (|Io|). We assumed that for 
a gluon probe [fl3[ we have the same relation between the gluon structure function and the measured cross 
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section. It should be stressed that Eq. (|18|) coincides with the definition of the gluon structure function in 
the DGLAP limit. 

Based on the model developed for the initial condition we will attempt to solve the evolution equation ( 
see Eq. (^)) and so provide reliable estimates for the structure functions in the LHC kinematic region. 



XIV. COMPARISON WITH OTHER MODELS: 



the models that are on the market fall into two classes: (i) models that introduce the saturation scale 
in addition to the separation one ( see our model and Refs. (IfjJi^^Tj ), and (ii) models that only have 
a separation scale (see Refs. |TJ,^ : ^G|,^8|). We discuss here only models that reproduce the same data or 
almost the same data as ours with a chi squared of approximately the same value. There are two such 
models^] 

Our model and Golec-Biernat and Wusthoff model jlo| are in the first class: both of them introduce a 
saturation scale and we will discuss the difference between them a little later. The Forshaw, Kelley and 



Shaw model |4q] suggested a different approach based on ideas of Donnachie and Landshoff model 1 14 3 



They use Eq. (Bj), and for the dipole cross section they suggest the sum of "soft" and "hard" contribution: 

<J dtpo ie(ri,W) = fs(r±)W Xs + f H {r± (r 2 ± W 2 ) Xh , (19) 



with the following limits of 

fs(r±) — > Const at r ± > 0.8-=- 1 fm ; (20) 
fH(r±) — ► r\ at r x < 028^- 0. 3 fm. (21) 

Therefore, in Eq. (M) one can recognize the "soft" Pomeron in the first term, while the second one at short 
distances ( r± <C r± ) depends on x as it should for the "hard" contribution. 

The success of this model shows us that the HERA data alone cannot distinguish between a rather small 
separation scale and saturation of the gluon density. In this model "hard" and "soft" are compatible at 
r± ~ 0.4 fm which is rather small. 

There is no difference in principle between the Golec- Wusthoff model and ours since both of them use the 
same physical picture that follows from the high density QCD. However, in our model the saturation scale 
can be calculated (see Eq. (||) ) and all our results have the correct DGLAP limit at short distances. On 
the other hand, Golec- Wusthoff model introduces the saturation scale as a phenomenological parameter, has 
an extremely simple form and reproduces the HERA data on the energy behaviour of ratio for crfg{ /o'tot- 
In other words, Golec- Wusthoff model investigates the extreme limit: all physics observed at HERA is 
related to a saturation scale. It is very encouraging that such a rudimentary idea does not contradict all 
experimental data from HERA. 

In Fig. [l8| we plot the ratio of our dipole cross section to Golec- Wusthoff one. More work is needed to 
understand how much of the disagreement stems from attractive features of our models such as correct 
matching with the DGLAP evolution equation, and how much of them can be explained by the fact that 
our model underestimates the value of SC at rather low Q 2 ss 1 2 GeV 2 . 



2 Model of Ref. [^7j also introduces a saturation scale but in a quite different way than in our model and in Golec- 
Biernat and Wusthoff model ^(J . However, this model is only in a very premature stage of development and has not 
described the HERA experimental data with a good chi squared. One can find the comparison of this model with 
other competing models in Refs. p^ , ^(i|j48| . 
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Fig. 18-a Fig.l8-b 
FIG. 18. The ratio of our dipole cross section to Golec-Biernat and Wusthoff dipole cross section at x = 10~ 3 
versus r\(GeV~ 2 } for total cross section ( Fig. 18-a) and for diffractive cross section (Fig. 18-b). 

It is incorrect to say that it is not possibile to distinguish one model from another. For example, for 
diffractive production it is important to take into account the diffractive dissociation in qqG channel. This 
channel can be considered as the production of a dipole but with a charge which is twice larger ( for N c > 1 
) than for a qq dipole. In spite of this fact in the two models ^,^8| the dipole for qqG - channel was treated 
in the same way as qq dipole ( in MFGS model Ji?]] it has not yet been considered ). In our model which 
is based on the correct matching with the pQCD result, the propagation of qqG - dipole is quite different 
than for a quark-anti quark pair. This means that these models should have different predictions for jet 
production in the diffractive dissociation processes. 



XV. OUR ANSWER: 



In our opinion HERA has reached a high density QCD domain because: (i) HERA data display a large 
value of the gluon structure function; (ii) no contradictions with the asymptotic predictions of high density 
QCD have been observed; and (iii) the numerical estimates of our model give a natural description of the 
size of the deviation from the usual DGLAP explanation. However, we know the main shortcoming of our 
statement, which is that most of the observed indications of high density effects can be described alternatively 
without assuming a saturation scale, but by matching "soft" and "hard" processes at rather short distances. 
We want to stress three aspects: 

1. The alternative approaches cannot describe or, perhaps, have not described all data and, in particu- 
lar, the energy dependence of the 3/^> photon and DIS production p5| . We found that a simultaneous 
analysis of this reaction and of the F2 slope cannot be done without employing shadowing corrections; 

2. The gluon saturation is not an additional postulate of pQCD, but follows from the QCD evolution 
equations in the high parton density kinematic region, and our model describes the solution to these 
equations; 

3. The gluon saturation leads to the simplest Golec-Biernat and Wusthoff model which gives an im- 
pressive description of the data. 

The fact that one can describe the HERA data without assuming gluon saturation is not very surprising 
as the saturation scale in the HERA kinematic region is Q 2 s {x) — 1 3GeV 2 while the typical scale for the 
"soft" Pomeron is also not very small and can be as large as 4GeV 2 |49[ . 

We hope that the simultaneous analysis of all HERA data , which we plan to undertake, will be consistent 
only with the saturation approach. In doing so we have to rely on the hdQCD evolution equations rather 
then on models. 
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XVI. PREDICTIONS FOR THERA: 



As we have pointed out it is rather dangerous to discuss predictions for lower x in our model, since the 
model should be replaced by the solution of the correct evolution equation of Eq. (||). However, in this 
section we summarize some of our predictions for as low x as x — > 10~ 6 since this region of x is under 
detailed discussion due to plans of the THERA extension of TESLA project |50[| . We would like to recall 
that our predictions underestimate the solution to the nonlinear evolution equation, but we can use them as 
a first estimate in our search for collective phenomena in high parton density system. 



The unitarity boundary and THERA kinematic region: 



We start from the prediction which, in principle, does not depend on the exact form of the correct evolution 
equation, namely, from the unitarity boundary for the F2 and xG(x,Q 2 ) structure functions. This problem 
has been considered in Ref. and it turns out that we have the following unitarity boundaries for F2 and 
xG(x,Q 2 ): 

dF 2 {x,Q 2 ) 1 2 2 . 

d\nQ* 3^2 ^ ' (22) 

d 2 xG(x,Q 2 ) 2 



< —oQR, (23) 



SlnQ 2 <91n(l/a:) tt 

where R is the nucleon radius. Actually, in Eq. ( p2| ) and Eq. ( p3| ) we need to R 2 —< b\ >p 2 and R 2 =< 
b 2 >xG, respectively (see Fig. 7). As it shown in Fig. 7, the interaction radii depend on energy and Q 2 but 
the dependence is so weak that we can neglect it in the first approximation. 



In Fig. |19| we plot Eq. (g3|) and the predictions for F 2 -slopes in GRV |32j,in MRS |3l| and in CTEQ |0| 
parameterizations. One can see that the i^-slope in the DGLAP evolution equations reaches the unitarity 
boundary at Q 2 < 3 GeV 2 for x = 1(T 5 and even for Q 2 < 10 GeV 2 for x = 1(T 6 . 

Fig. ^0| shows the comparison of the DGLAP evolution equations with the unitarity boundary for the case 
of gluon structure function. In Fig. Eq. ( p3| ) is modified using the DGLAP evolution equation in the 
region of low x, namely, 

d 2 xG DGLAP (x,Q 2 ) N c a s (Q 2 ) ^dglap, q2) (9A] 

ainQ 2 dHl/x) - " xG > ■ ^ 

Using Eq. ( p^ ) one can rewrite Eq. J23| ) in the form 

xG DGLAP (x,Q 2 ) < 2 Q 2 R 2 . (25) 

7r N c a s (Q z ) 
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CTEQ5 

unitarity boundary 



FIG. 19. The unitarity boundary for F^-slope and predictions in GRV,MRS and CTEQ parameterizations for the 
solution of the DGLAP evolution equations in THERA kinematic region. 

Based on Fig. ^0]we are not optimistic of finding gluon saturation for the unitarity constrains on the gluon 
structure function in THERA kinematic region. The -FVslope at lower x (x ss 10 -6 ) will give more reliable 
limit on the gluon saturation. 
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FIG. 20. The unitarity boundary for xG(x, Q 2 ) and predictions in GRV,MRS and CTEQ parameterizations of the 
solution of the DGLAP evolutions equations in THERA kinematic region. 

However, the calculations, given in Fig. 6 using both our model and the solution to the Eq. (|J), suggest 
that the corrections to at least xG is rather large in the THERA kinematic region. 



dF 2 
ainQ 2 

The THERA kinematic region provides a new opportunity to check the value of SC using the F2 slope. At 
least, the THERA data would help us to differentiate between two approaches: our approach, based on the 
gluon saturation at low x and the approach based on the matching between "soft" and "hard" interaction 
with a sufficiently large momentum scale for "soft" contribution. Fig. |2l] shows the predictions for the 
F 2 -slope in our approach ( GLMN in figure ) and the Donnachie-Landshoff approach pffl. The difference is 
rather large and can be measured. 
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FIG. 21. The prediction for the Fi-slope in our approach (GLMN) and in Donnachie-Landshoff approach (DL). 



J/SE' production: 



Fig. 22 shows that the difference between our predictions and the DGLAP increase at the THERA 
kinematic region, which leads us to speculate that we will see gluon saturation at THERA. 



Q 2 =0 



Q 2 = 10. 1 GeV 2 




-3 200 




FIG. 22. Predictions for J/fy production in THERA kinematic region together with available experimental data. 
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Maxima in the ratios and the higher twist contributions 



Fig. [I?] shows that we expect a characteristic behaviour of the ratio /F® as a function of Q 2 . It turns 
out to be wider than at higher x and the maximum occurs at the value of Q^ax ~ 6 7 GeV 2 . Such a large 
value of Q 2 nax makes the calculations reliable and we expect that the measurement of this observable in the 
THERA kinematic region will help us to extract the value of the saturation scale from the experimental 
data. In THERA kinematic region ( at x « 10~ 5 ) we expect that the higher twist contributions to be of the 
same order as the leading twist ones at sufficiently high value of Q 2 , namely, for Fl and for F® this value 
of Q 2 is about 5 ~ 7 GeV 2 while for F® it is even larger, about 15 GeV 2 . The fact that the higher twist 
contributions become visible at large value of Q 2 lead to a visible violation of the usual DGLAP evolution 
approach which should be possible to measure. 



ACKNOWLEDGMENTS 



The authors are very much indebted to our coauthors and friends with whom we discussed our approach 
on a everyday basis Ian BalitskyJochen Bartels , Krystoff Golec Biernat, Larry McLerran, Dima Kharzeev, 
Yuri Kovchegov and Al Mueller for their help and fruitful discussions on the subject. E.G. , E. L. and U.M. 
thank BNL Nuclear Theory Group and DESY Theory group for their hospitality and creative atmosphere 
during several stages of this work. 

This paper is the fruit of the useful discussions during Amirim meeting and it is a pleasure for us to 
thank all participants of the Amirim workshop: H. Abramowicz, J. Bartels, L. Frankfurt, K. Golec-Biernat, 
E. Gurvich, H. Jung, S. Kananov, H. Kowalsky, A. Kreisel, A. Levy, P. Marage, M. McDermott, P. Saull, 
S. Schlenstedt, G. Shaw, M. Strikman, J. Whitmore . 

This research was supported in part by the BSF grant # 9800276 and by Israeli Science Foundation, 
founded by the Israeli Academy of Science and Humanities. 



[1] L.V. Gribov, E.M. Levin and M.G. Ryskin, Phys. Rep 100 (1983) 1; 
A.H. Mueller and J. Qiu, Nucl. Phys. B268 (1986) 427; 

L. McLerran and R. Venugopalan,P/iys. Rev. D49 (1994) 2233,3352, 50 (1994) 2225, 53 (1996) 458, 59 (1999) 
094002. 

[2] V.N. Gribov, Sov. Phys. JETP 30 (1970) 709. 

[3] E. Levin and M.G. Ryskin, Phys. Rep. 189 (19267) 1990; 

J.C.Collins and J. Kwiecinski, Nucl. Phys. B335 (1990) 89; 

J. Bartels, J. Blumlein and G. Shuler, Z. Phys. C50 (1991) 91; 

E. Laenen and E. Levin, Ann. Rev. Nucl. Part. Set. 44 (1994) 199 and references therein; 
AL. Ayala, M.B. Gay Ducati and E.M. Levin, Nucl. Phys. B493 (1997) 305, B510 (1990) 355; 
la. Balitsky, Nucl. Phys. B463 (1996) 99; 

Yu. Kovchegov, Phys. Rev. D54 (191996) 5463, D55(1997) 5445, D60(2000) 034008, D61 (2000)074018; 
A.H. Mueller, Nucl. Phys. B572(2000)227, B558 (1999) 285; 
Yu. V. Kovchegov, A.H. Mueller, Nucl. Phys. B529 (1998) 451; 
E. Levin and K. Tuchin, Nucl. Phys. B573(2000) 833. 
[4] A.H. Mueller, Nucl. Phys. B415 (1994) 373. 

[5] J. Jalilian-Marian, A. Kovner, L. McLerran and H. Weigert, Phys. Rev. DD55 (1997) 5414; 
J. Jalilian-Marian, A. Kovner and H. Weigert, Phys. Rev. D59 (1999) 014015; 
J. Jalilian-Marian, A. Kovner and H. Weigert, Phys. Rev. D59 (1999) 014015; 

J. Jalilian-Marian, A. Kovner, A. Leonidov and H. Weigert, Phys. Rev. D59 (1999) 014014,034007, Erratum- 
ibid. Phys. Rev. D59 (1999) 099903; 

A. Kovner, J.Guilherme Milhano an d H. Weigert, OU TP-00-10P,NORDITA-2000-14-HE, friep-ph/0004014| ; 
H. Weigert, NORDITA-2000-34-HE, |hep-ph/0004044|. 



27 



Ia. Balitsky, Nud.Phys. B463 (1996) 99; Yu. Kovchegov, Phys. Rev. D60 (191000) 034008. 

A.L. Ayala, M.B. Gay Ducati and E.M. Levin, Nucl. Phys. B493 (1997) 305, B510 (1998) 355. 

L.McLerran, " Renormalization Group Equations, Saturation, and the Colored Glass Condensate", talk given at 

WS "QCD in a nuclear environment", August 3-5, Regensburg, Germany. 

Yu. V. Kovchegov, Phys. Rev. D61 (2000) 074018; 

E. Levin and K. Tuchin, Nucl.Phys. B573 (2000) 833; 

M. Braun, Eur. Phys. J. C16 (2000) 337. 

E. Levin and L. Frankfurt, JETP Lett. 2 (1965) 65; H. J. Lipkin and F. Scheck, Phys. Rev. Lett. 16 (1966) 71. 
A. Zamolodchikov, B. Kopeliovich and L. Lapidus, JETP Lett. 33 (1981) 595; 
E.M. Levin and M.G. Ryskin, Sov. J. Nucl. Phys. 45 (1987) 150. 
A. H. Mueller, Nucl. Phys. B335 (1990) 115. 

A. Donnachie, P. V. Landshoff,P%i. Lett. B296 (1992) 227; B437(1998) 408 a nd references the rein. 



G. Cvetic, D. Schildknecht and A. Shoshi.Acta Phys.Polon. B30 (1999) 3265, frep-ph/9910379| and references 
therein. 

E. Shuryak and T. Schafer, Ann. Rev. Nucl. Part. Set. 47 (1997) 359 and references therein. 
E. Gotsman, E. Levin, U. Maor and E. Naftali, Eur. Phys. J. CIO (1999) 689; 
E. Gotsman, E. Levin and U. Maor, Eur. Phys. J. C5 (1998) 303. 
N.N. Nikolaev and B.G. Zakharov, Z. Phys. C49 (1991) 607; 

E.M. Levin, A.D. Martin, M.G. Ryskin and T. Teubner, Z. Phys. C74 (1997) 671. 

E. Gotsman, E. Levin and U. Maor, Phys. Lett. B379 (1990) 186. 

E. Gotsman, E. Levin, U. Maor and E. Naftali, Nucl.Phys. B539 (1999) 535; 

E. Gotsman, E. Levin and U. Maor, Phys. Lett. B425 (1998) 369. 

E. Gotsman, E. Levin, U. Maor and E. Naftali, Eur. Phys. J. C14 (2000) 511. 

E. Gotsman, E. Levin and U. Maor, Nucl.Phys. B493 (1997) 354. 

DD 

E. Gotsman, E. Levin, M. Lublinsky, U. Maor and K. Tuch in, " Energy dependence of — — in DIS and shadowing 



corrections", TAUP-2639-2000, Jul 2000, tiep-ph/0007261 



E. Gotsman, E. Levin and U. Maor, Phys. Lett. B403 (1997) 120; Nucl.Phys. B464 (1996) 251. 

E. Gotsman, E. Ferreira, E. Levin, U. Maor and E. Naftali, " Screening corrections in DIS at low Q 2 and x" , 

Talk given at 30th Internat ional Conference on High-Energy Physics (ICHEP 2000), Osaka, Japan, 27 Jul - 2 



Aug 2000, hep-ph/0007274 



V.A. Abramovsky, V.N. Gribov and 0,V. Kancheli, Sov. J. Nucl. Phys. 18 (1974) 308. 

J. Bartels, Phys. Lett. B298 (1993) 204, Z. Phys. C60 (1993) 471; 

E.M. Levin, M.G. Ryskin and A.G. Shuvaev, Nucl. Phys. B387 (1992) 589. 

E. Gotsman, E. Levin, U. Maor, L. McLerran and K. Tuchin, "Higher twists an d maxima for DIS on nuclei in 
high density QCZ)",TAUP-2638-200, BNL-NT-00-19, Jul 2000, friep-ph/0007256 



E. Gotsman, E. Levin, M. Lublinsky, U. Maor and K. Tuchin, "Shadowing corrections and diffractive production 



in DIS on nuclei", TAUP-2605-99, Nov 1999, [hep-ph/9911270 . 

Yu. V. Kovchegov and L. McLerran, Phys. Rev. D60 (1999) 054025. 

H. L. Lai et al., Eur. Phys. J. C12 (2000) 37; 

A. D. Martin, R. G. Roberts, W. J. Stirling and R. S. Thorne, Eur. Phys. J. C4 (1998) 463. 
M. Gluck, E. Reya and A. Vogt, Eur. Phys. J. C5 (1998) 461. 

HI collaboration: S. Aid et al, Nucl. Phys. B470 (1996) 3; C, Adolf et al., Nucl. Phys. B497 (1997) 3; 
ZEUS collaboration : M. Derrick et a l., Z. Phys. C72 (1996) 399; J. Breitweg et al., Phys. Lett. B407 (1997) 



432; DESY-00-071, tiep-ex/000501£ 



A. M. Cooper-Sarkar, R. C. E. Devenish and A. De Roeck, Int. J. Mod. Phys. A13 (1998) 33; 
H. Abramowicz and A. Caldwell, Rev. Mod. Phys. 71 (1999) 1275. 

M. Klein (HI collaboration), "Structure Fun ctions in Deep Ine lastic Lepton-Nucleon Scattering" Talk at Lepton- 



Photon Symposium, Stanford, August 1999, hep-ex/0001059 



B. Foster (ZEUS Collaboration), Invited talk, Royal Society Meeting, London, May 2000; 
ZEUS collaboration: J. Breitweg et al., Phys. Lett. B487 (2000) 53. 

A. Donnachie and P.V. Landshoff: Phys. Lett. B437 (1998) 408; Phys. Lett. B470 (1999) 243. 
HI Collaboration: C. Adloff et al., Z. Phys. C76 (1997) 613; 
ZEUS Collaboration: J. Breitweq et al., Eur. Phys. J. C6 (1999) 43. 
E. Levin and M. Wusthoff, Phys. Rev. D50 (1994) 4306; 
M. Wusthoff and A.D. Martin, J. Phys. G2 (1999) R309. 

K. Golec-Biernat and M. Wusthoff, Phys. Rev. D59 (1999) 014017; D60 (1999) 114023; K.Golec-Biernat,Talk 
at 8th International Workshop on Deep Inelastic Scattering and QCD (DIS 2000), Liverpool, England, 25-30 



28 



[41] 
[42] 

[43] 

[44] 
[45] 



[46] 
[47] 

[48] 
[49] 



[50] 
[51] 



Apr 2000 |iep-ph/000608c| . 

ZEUS collaboration: J. Breitweg et al, Eur. Phys. J. C14 (2000) 213. 

HI Collaboration: C. Adloff et &\.,Phys. Lett. B483 (2000) 36. 

J. Bartels, A. De Roeck and H. Lotter, Phys. Lett. B389 (1996) 742; 

J. Bartels, A. De Roeck, C. Ewerz and H. Lo tter, "The gamm a* gamma* total cross section and the BFKL 

pomeron at the 500-GeV e+ e- linear collider", tiep-ph/971050C ; 

S. J. Brodsky, F. Hautmann and D. E. Soper, Phys. Rev. Lett. 78 (1997) 803. 

PLUTO collaboration: C. Berger et al., Phys. Lett. B149 (191984) 42, Z. Phys. C26 (1984) 353; 

PC/Two Gamma Collaboration: H. Aihara et al., Phys. Rev. D41 (1990) 2667;. D. Bintingeret al., Phys. Rev. 

Lett. 54 (1985) 763; 

OPAL Collaboration: G. Abbiendi et al., Eur. Phys. J. C14 (2000) 1; 
L3 Collaboration: M. Acciarri et al., Phys. Lett. B408 (1997) 450; 

E. Gotsman, E. Levin, L. McLerran and K. Tuchin, " Higher twists and maxima for DIS on proton in high density 

QCD", TAUP-2644-2000, BNL-NT-00-20, Aug 2000, [hep-ph/0008280| . 

B. Badelek and J. Kwiecinski, Phys. Lett. B418 (1998) 229 and references therein; 

J. Bartels, K. Golec-Biernat and K. Peters, " An estimate of higher twist at small xb and low Q 2 based upon a 
saturation model", DESY-00-038, Mar 2000, [hep-ph/0003042| . 

M.F. McDermott, " The dipole picture of small x physics ( a summary of the Amirim meeting)", hep-ph/0008260 
M. McDermott, L. Frankfurt, V. Guzey and M. Strikman, "Unitarity and the QCD-improved dipole picture", 



hep-ph/9912547 



J. R. Forshaw, G. Kerley and G. Shaw, Phys. Re. D60 (1999) 074012, Nucl. Phys. A675 (2000) 80. 
D. Kharzeev and E. Levin, Nucl.Phys. B578 (2000) 351; 

Yu. Kovchegov, D. Kharzeev and E. Levin, " QCD instantons and the soft Pomeron ",BNL-NT-00-18, TAUP- 
2637-2000, Jul 2000, |hep-ph/0007182| ; 

B.Z. Kopeliovich, I.K. Potashnikova, B. Povh and E. Predazzi, Phys. Rev. Lett. 85 (2000) 507. 

Max Klein, "THERA-electron proton scattering at ^/s ~ ITeV", talk, given at DIS'2000, Liverpool, April 

25-30,200; and web-page: www.ifh.de/thera/. 

A.L. Ayala Filho, M.B. Gay Ducati and E.M. Levin,Phys.Lett. B388 (1996) 188. 



29 



